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ABSTRACT 

To investigate the black hole mass vs. stellar velocity dispersion (Mbh — &*) relation of active 
galaxies, we measured the velocity dispersions of a sample of local Seyfert 1 galaxies, for which we 
have recently determined black hole masses using reverberation mapping. For most objects, stellar 
velocity dispersions were measured from high signal-to-noise ratio optical spectra centered on the Ca II 
triplet region (~ 8500 A), obtained at the Keck, Palomar, and Lick Observatories. For two objects, 
in which the Ca II triplet region was contaminated by nuclear emission, the measurement was based 
on high-quality _ff-band spectra obtained with the OH-Suppressing Infrared Imaging Spectrograph at 
the Keck-II Telescope. Combining our new measurements with data from the literature, we assemble 
a sample of 24 active galaxies with stellar velocity dispersions and reverberation-based black hole 
mass measurements in the range of black hole mass 10 6 < Mbh/M© < 10 9 . We use this sample to 
obtain reverberation mapping constraints on the slope and intrinsic scatter of the Mbh — relation 
of active galaxies. Assuming a constant virial coefficient / for the reverberation mapping black hole 
masses, we find a slope (3 = 3.55 ± 0.60 and the intrinsic scatter a lnt = 0.43 ± 0.08 dex in the relation 
log(MBH/M ) = a + ;31og(cT*/200 km s~ 4 ), which are consistent with those found for quiescent 
galaxies. We derive an updated value of the virial coefficient / by finding the value which places the 
reverberation masses in best agreement with the Mbh — c* relation of quiescent galaxies; using the 
quiescent A/bh — relation determined by Giiltekin et al. we find log / = 0.72lg'°g with an intrinsic 
scatter of 0.44±0.07 dex. No strong correlations between / and parameters connected to the physics of 
accretion (such as the Eddington ratio or line-shape measurements) are found. The uncertainty of the 
virial coefficient remains one of the main sources of the uncertainty in black hole mass determinations 
using reverberation mapping, and therefore also in single-epoch spectroscopic estimates of black hole 
masses in active galaxies. 

Subject headings: galaxies: active — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: 
Seyfert 



1. INTRODUCTION 

In the past decade, supermassive black holes have been 
found to be virtually ubiquitous at the centers of local 
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quiescent galaxies (for reviews, see Kormendy 2004; Fer- 
rarese & Ford 2005). Remarkably, the mass of the central 
black hole (Mbh) correlates with the global properties of 
their host spheroids, primarily the stellar velocity disper- 
sion (cr*; Ferrarese et al. 2000; Gebhardt et al. 2000a). 
This Mbh — c* relation is believed to be of fundamen- 
tal importance, providing one of the strongest empirical 
links between galaxy formation/evolution and nuclear ac- 
tivity. Measuring the slope and scatter of the relation, 
and any possible dependence on additional parameters, 
is essential to make progress in our understanding of the 
co-evolution of galaxies and black holes. 

One key open question is whether active galactic nu- 
clei (AGNs) share the same A/bh — c* relation as quies- 
cent galaxies, as would be expected if the relation were 
universal and nuclear activity were just a random tran- 
sient phase. Initial studies using several reverberation- 
mapped Seyfert 1 galaxies show that the Mbh — c* re- 
lation of active galaxies is consistent with that of qui- 
escent galaxies (Gebhardt et al. 2000b; Ferrarese et al. 
2001). Other studies using larger samples with single- 
epoch black hole masses suggest that the slope of the 
relation may be different than that of quiescent galaxies, 
albeit at low significance (e.g., Greene & Ho 2006; Shen 
et al. 2008). However, these later studies rely on an indi- 
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rect estimator of black hole mass, based on the kinemat- 
ics and size of the broad-line region (BLR) as inferred 
from single-epoch spectra. The method works as follows: 
the broad-line profile gives a characteristic velocity scale 
AV (measured from either line dispersion, <7n nc , or the 
full- width at half-maximum intensity, Vfwhm); the av- 
erage size of the BLR (-Rblr) is determined from the 
empirical correlation with optical luminosity (Wandel et 
al. 1999; Kaspi et al. 2005; Bentz et al. 2009a); and the 
black hole mass is obtained as 



M BH = / 



(AU) 2 i? B LR 

G 



(1) 



where G is the gravitational constant and / is a virial 
coefficient that depends on the kinematics and the geom- 
etry of the BLR. Although single-epoch mass estimates 
are believed to be accurate to ~ 0.5 dex (e.g., Vester- 
gaard & Peterson 2006), it is hard to quantify precisely 
their uncertainty and any possible biases (Marconi et al. 
2008; Onken et al. 2009). Thus, the slope and intrinsic 
scatter of the Mbh — cr* relation of active galaxies remain 
highly uncertain. 

Multi-epoch data provide direct measurements of the 
BLR size via reverberation mapping time lags (-Rblr = 
cr, where c is the speed of light and r is the measured re- 
verberation time scale), and more secure measurements 
of the broad- line kinematics as determined from the vari- 
able component of the broad lines (Peterson et al. 2004). 
Onken et al. (2004) combined reverberation black hole 
masses with measurements of stellar velocity dispersion 
for 14 objects to measure the slope of the Mbh — cr* re- 
lation, showing that AGNs and quiescent galaxies lie on 
the same A/bh — cr* relation for an appropriate choice 
of the virial coefficient (/ = 5.5 ± 1.8)0 Remarkably, 
the scatter in A:Tbh relative to the Mbh — cr* relation is 
found to be less than a factor of 3. Unfortunately, their 
sample is very small (14 objects) and covers a limited dy- 
namic range (6.2 < logMBH < 8.4), making it difficult 
to simultaneously measure the intrinsic scatter and the 
slope and investigate any trends with black hole mass, or 
with other properties of the nucleus or the host galaxy 
(cf. Collin et al. 2006). 

Under the assumption that the relations should be the 
same for active and quiescent galaxies, the Onken et al. 
(2004) study provides some information on the geome- 
try of the BLR, an absolute normalization of reverbera- 
tion black hole mass, and an upper limit to the intrinsic 
scatter in the virial coefficient. Since the reverberation- 
mapped AGN sample is the "gold standard" used to cal- 
ibrate all single-epoch mass estimates (e.g., Woo & Urry 
2002; Vestergaard et al. 2002; McLure & Jarvis 2002; 
Vestergaard & Peterson 2006; McGill et al. 2008; Shen 
et al. 2008), this comparison is effectively a crucial link 
in establishing black hole masses for all AGNs across the 
universe and for all evolutionary studies (e.g., McLure & 
Dunlop 2004; Woo et al. 2006, 2008; Peng et al. 2006a,b; 
Netzer et al. 2007; Bennert et al. 2010; Jahnke et al. 
2009; Merloni et al. 2010). 

In this paper, we present new measurements of host- 

15 Note that the virial coefficient, (/), is 5.5 when the line disper- 
sion (oi; ne ) is used. If Vfwhm i s used, (/} has to be properly scaled 
depending on the Vfwhm /""line ratio. See Onken et al. (2004) and 
Collin et al. (2006) for details. 



galaxy velocity dispersion obtained from deep, high- 
resolution spectroscopy at the Keck, Palomar, and Lick 
Observatories. The new measurements are combined 
with existing data (e.g., Onken et al. 2004; Nelson et al. 
2004; Watson et al. 2008) and with recently determined 
black hole masses from reverberation mapping (Bentz et 
al. 2009b) to construct the Mbh ~ cr* relation of broad- 
lined AGNs for an enlarged sample of 24 objects. For 
the first time, we are able to determine simultaneously 
the intrinsic scatter and the slope of the Mbh — cr* rela- 
tion of active galaxies. Also, by forcing the slope of the 
Mbh - cr* relation to match that of quiescent galaxies, we 
determine the average virial coefficient and its nonzero 
intrinsic scatter. Finally, we study the residuals from 
the determined Mbh — cr* relation to investigate possible 
trends in virial coefficient with properties of the active 
nucleus. 

The paper is organized as follows. In §2, we describe 
the observations and data reduction of the optical and 
the near-infrared (IR) data. New velocity dispersion 
measurements are presented in §3, along with a sum- 
mary of previous measurements in the literature. In §4, 
we investigate the Mbh — cr* relation of the present-day 
AGNs, and determine the virial coefficient in measuring 
Mbh- Discussion and summary follow in §5 and §6, re- 
spectively. 

2. OBSERVATIONS 

The Lick AGN Monitoring Project (LAMP) was de- 
signed to determine the reverberation time scales of a 
sample of 13 local Seyfert 1 galaxies, particularly with 
low black hole masses (< 1O 7 M ). NGC 5548, the best- 
studied reverberation target, was included in our Lick 
monitoring campaign to test the consistency of our re- 
sults with previous measurements. The 64-night spectro- 
scopic monitoring campaign, along with nightly photo- 
metric monitoring, was carried out using the 3-m Shane 
reflector at Lick Observatory and other, smaller tele- 
scopes. The detailed photometric and spectroscopic re- 
sults are described by Walsh et al. (2009) and Bentz et 
al. (2009b, 2010), respectively. In summary, nine objects 
including NGC 5548 showed enough variability in the op- 
tical continuum and the H/3 line to obtain the reverber- 
ation time scales. This significantly increases the size of 
the reverberation sample, particularly for this low-mass 
range. 

For the LAMP sample of 13 objects, we carried out 
spectroscopy using various telescopes, to measure stellar 
velocity dispersions. Here, we describe each observation 
and the data-reduction procedures. The date, total ex- 
posure time, and other parameters for each observation 
are listed in Table [TJ 

2.1. Optical Data 
2.1.1. Palomar Observations 

Observations of six galaxie s were obtained w ith the 
Double Spectrograph fDBSP: lOke fc Gundfl98l at the 
Palomar Hale 5-m telescope (P200). A 2"-wide slit was 
used, along with the D68 dichroic. On the red side of the 
spectrograph, we used a 1200 lines mm -1 grating blazed 
at 9400 A, covering the wavelength range 8330-8960 A 
at an instrumental dispersion of crj ~ 30.4 km s _1 . 

Each galaxy was observed with the slit oriented ap- 
proximately at the parallactic angle (Filippenko 1982) for 
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TABLE 1 
Observation Log 



Galaxy 


- — 

Telescope / 


UT Date 


Exposure Time 


PA 


Airmass 




S/N 




Instrument 




0) 


(deg) 






Arp 151 


P200/DBSP 


2003-01-28 


5400 


210 


1.10 


79 


IC 1198 (Mrk 871) 


Keck/ESI 


2008-03-02 


900 


285 


1.25 


87 


IC 4218 


P200/DBSP 


2003-06-02 


3600 


5 


1.23 


61 


MCG-06-30-15 


P200/DBSP 


2003-06-01 


5400 


10 


2.70 


112 


Mrk 1310 


Keck/ESI 


2008-03-02 


1800 


323 


1.19 


84 


Mrk 142 


Keck/OSIRIS 


2009-05-03 


3600 





1.19 


17 


Mrk 202 


Keck/ESI 


2008-03-02 


1200 


204 


1.32 


73 


Mrk 290 


Lick/Kast 


2008-04-13 


3600 


168 


1.07 


27 


NGC 4253 (Mrk 766) 


P200/DBSP 


2001-06-26 


1200 


220 


1.60 


63 




Keck/OSIRIS 


2009-05-03 


3600 





1.03 


67 


NGC 4748 


Keck/ESI 


2004-02-17 


1200 


25 


1.23 


160 




Keck/OSIRIS 


2009-05-04 


5400 





1.37 


61 


NGC 5548 


P200/DBSP 


2003-06-01 


1800 


59 


1.06 


100 


NGC 6814 


P200/DBSP 


2003-06-01 


3600 





1.38 


165 


SBS 1116+583A 


Keck/ESI 


2008-03-02 


1800 


216 


1.40 


53 



Note. — For OSIRIS observations, the position angle (PA) refers to the direction of the long 
axis of the IFU; for all other observations it refers to the slit PA. The S/N refers to the signal- 
to-noise ratio per pixel in the extracted spectra, at ^8400-8700 A for the Ca II triplet spectral 
region, or at ~1. 47-1. 61 /an for the ff-band spectra. 



the midpoint of the observation. Typically two or three 
exposures were taken for each galaxy to aid in cosmic-ray 
removal. Flux standards and a range of velocity tem- 
plate stars (primarily K-type giants) were observed dur- 
ing each night. 

2.1.2. Lick Observations 

Mrk 290 was observed with the Kast Double Spectro- 
graph (Miller & Stone 1993) at the Shane 3-m telescope 
at Lick Observatory on 2008 April 13 (UT dates are 
used throughout this paper), during our AGN monitor- 
ing campaign. For these observations, we used the D55 
dichroic and the 830/8460 grating on the red side of the 
spectrograph, covering the wavelength range 7570-9620 
A at a scale of 1.7 A pixel" 1 . A 2"-wide slit was used 
and oriented at the parallactic angle for the midpoint of 
the exposure sequence, and the instrumental dispersion 
was at ~ 59 km s -1 at 8600 A. We obtained four 900 s 
exposures of Mrk 290 with this setup, along with short 
exposures of three velocity template stars and the flux 
standard star BD+28°4211. 

2.1.3. Keck Observations 

Five galaxies were observ ed with the Echelle te Spec- 
trograph and Imager (ESI; iSheinis et al.l l2002t ) at the 
Keck-II 10-m telescope. In ESI echellette mode, the ob- 
servations cover 3900-11000 A in 10 spectral orders at a 
scale of ~11.4 km s -1 pixel -1 . We used a O'.'75-wide slit, 
giving an instrumental dispersion of Oi ~ 22 km s . 

The slit was oriented at the parallactic angle in all 
observations. During twilight of each night, we observed 
flux standards and several velocity template stars with 
spectral types ranging from G8III to K5III. 

2.1.4. Reductions 

We reduced the DBSP and Kast data using a series 
of IRAFr^l scripts. The standard spectroscopic data- 
reduction procedure, including bias subtraction, flat 

16 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 



fielding, wavelength calibration, spectral extraction, and 
flux calibration was perform ed for each data set. Opti- 
mal extraction (Home 1986) was used for obtaining one- 
dimensional spectra to achieve maximal signal-to-noise 
ratio (S/N) on the stellar features. The typical extrac- 
tion radius was ~ 2 — 3", which corresponds to a physical 
radius of ~ 1 kpc. 

The ESI observations require some special prepara- 
tion and calibration before they can be combined and 
used for kinematic analysis. The calibration process 
undertaken includes the following steps: bias subtrac- 
tion, flat-fielding, cosmic-ray rejection, wavelength cali- 
bration, rectification, and sky subtraction. These steps 
were performed by the IRAF package EASI2D, which 
was developed by David J. Sand and Tommaso Treu 
(jSand et al.l 120041 ) for easy extraction of echelle orders. 
Approximate flux calibration was performed using ESI 
response curves measured during photometric nights. 
During wavelength calibration the spectra were rebinned 
to uniform steps in log(A) corresponding to 11.4 km s , 
which is close to the native pixel scale of ESI, thus min- 
imizing covariance between pixels and loss of resolution 
due to rebinning. The root-mean square (rms) residu- 
als in the wavelength solution are much smaller than a 
pixel (typically 5%) and therefore negligible. Spectral 
resolution was measured from night-sky lines and from 
wavelength-calibration lines taken in the same configura- 
tion and reduced through the same pipeline to include all 
instrumental effects. As for the DBSP and Kast data, op- 
timal extraction was used for obtaining one-dimensional 
spectra. The typical extraction radius of the ESI spectra 
was ~ 1 — 1.5", which corresponds to a physical radius 
of ~ 0.5 kpc at the mean redshift of the observed AGNs. 

2.2. Near-Infrared Data 

2.2.1. Keck Observations 

For three galaxies, we obtained ii-band spectra us- 
ing the integral-field unit (IFU) OH-Suppressing Infrared 

for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Imaging Spectrograph (OSIRIS; lLarkin et al.l I2006D at 
the Keck-II telescope, operated with t he laser guide-star 
adap tive optics (LGS-AO) system (|Wizinowich et all 
2006). Using the high angular resolution in the AO- 
corrected data, our goal was to extract the light from 
an annulus around the nucleus if necessary, by excluding 
light from the cent ral, AGN-dom i nated region, similar to 
what was done bv lWatson et al.l (|2008t ). 

The OSIRIS observations were obtained on 2009 May 
3 and 4, during the first half of each night. We used 
the ff broadband (Hbb) hlter and the O'.'l lenslet scale, 
giving a field of view of 1'.'6 x 6'.'4, with spectra covering 
the range 1473-1803 nm at a spectral resolving power 
of R w 3800. The long spatial direction of the IFU was 
oriented north-south for all observations. 

The observing sequence for each galaxy consisted of 
sets of four exposures. In three exposures the nucleus 
of the galaxy was placed on the IFU. Then one sky ex- 
posure with an offset of 21" to the southeast was taken. 
The three on-source exposures in each set were dithered 
northward by 0'.'6 between exposures. Each individual 
exposure had a duration of 600 s for NGC 4748, or 300 s 
for Mrk 142 and Mrk 766, and total on-source exposure 
times for each source are listed in Table Q] Immediately 
preceding and following the observation sequence for each 
galaxy, we observed an A0V star for telluric correction 
and flux calibration. The A0V stars were observed using 
a sequence of two on-source and one off-source exposure, 
with typical exposure times of 5-10 s. In order to ob- 
tain unsaturated exposures of these bright stars, the AO 
loops were not closed during the star observations. On 
each night, velocity template stars were also observed, 
following the same observing sequence as used for the 
telluric correction stars. We observed several K and M- 
type giant stars as velocity templates, and A0V telluric 
correction stars were also observed close in time and air- 
mass to the velocity template stars. 

The nucleus of each galaxy was used as a tip-tilt refer- 
ence source (V ~ 15 — 16 mag) for the LGS-AO system, 
and the laser was propagated on- axis for each dithered 
exposure. It is difficult to quantify the image quality or 
determine the Strehl ratio in the galaxy exposures, since 
the galaxy nuclei are not point sources. However, the 
typical LGS-AO performance for tip-tilt point sources 
for these similar magnitudes yields if - band Strehl ratios 
of 40% to 50% (jvan Dam et al.ll2006l l20f)l . Therefore, 
given that the quality of both nights was excellent, we 
conservatively estimate the ff-band Strehl ratio to be 
~15%-20% for our galaxy observations. 

Prior to this OSIRIS run, the OSIRIS detector temper- 
ature had been rising from its previous operating tem- 
perature of 68 K. During these two nights, the detec- 
tor temperature was steady at 74 K, which resulted in 
a higher level (~40%) of dark current than normal for 
OSIRIS. The higher dark current did not significantly 
impact our observations, since our exposures were not 
dark or background limited. In the afternoon, multiple 
dark frames were taken with exposure times matching 
those of our observations, and were used as part of the 
sky-subtraction procedure. In addition, a new calibra- 
tion file for the observing mode with the Hbb filter and 
O'.'l scale was taken at a similar operating temperature 
to ensure a clean spectral extraction process. 



2.2.2. Reductions 

The data were reduced using the OSIRIS data- 
reduction pipeline software (jKrabbe et al.l [2004). In 
summary, the reduction steps included dark subtraction, 
cosmic-ray cleaning and glitch identification, extraction 
of the spectra into a data cube, wavelength calibration, 
sky subtraction using the offset sky exposures, correction 
for atmospheric dispersion, and telluric correction using 
the extracted spectra of the A0V stars. The OSIRIS 
pipeline includes an implementatio n of the scal ed sky- 
subtraction algorithm described by|Davics (2007?) , which 
provides a more accurate subtraction of the strong H- 
band airglow lines than a direct subtraction of the offset- 
sky exposure. However, the current implementation of 
the scaled sky routine left noticeable residuals around 
the wavelengths of the strongest blended OH features; 
we further discuss these sky residuals below. Following 
the sky subtraction and telluric correction, individual ex- 
posures were aligned and co-added to create a final data 
cube for each galaxy. 

We extracted 1-dimensional spectra from each data 
cube using rectangular extraction regions. While we 
had originally planned to extract an annular region 
around the nucleus in order to avoid the central, AGN- 
dominated region, we found that in the reduced data the 
-ff-band stellar absorption features were visible even in 
the central region of the galaxy, and it did not prove 
necessary to exclude the central region from the spectral 
extractions. We experimented with a variety of different 
extraction-region sizes for each galaxy. Since the galaxy- 
light profiles were very centrally concentrated, the ex- 
tracted spectra did not change significantly beyond some 
extraction aperture size, and in the end we used extrac- 
tion apertures that were large enough to contain the ma- 
jority of the detected galaxy light without including un- 
necessarily large amounts of noise from the surrounding 
regions. The extraction aperture sizes were 1'.'4 x 2'.'4, 
1'.'3 x 3'.'4, and 1'.'3 x 3'.'4, respectively, for Mrk 142, Mrk 
766, and NGC 4748. One-dimensional spectra were sim- 
ilarly extracted from the reduced data cubes for the ve- 
locity template stars. 

As a result of the imperfect sky subtraction mentioned 
previously, spectral regions around the wavelengths of 
the brightest OH emission blends were contaminated by 
noticeable residuals that appeared to affect all pixels on 
the IFU in a similar fashion. In order to correct for these 
additive residuals, for each galaxy we extracted spectra 
from rectangular regions at the outer edges of the IFU 
where there was very little galaxy continuum detected, 
with the same number of pixels as the nuclear extraction, 
and then subtracted these nearly blank-sky spectra from 
the nuclear spectra. This correction procedure effectively 
removed the night-sky emission residuals from the galaxy 
spectra. The procedure was applied for NGC 4748 and 
Mrk 142, but was not necessary for Mrk 766 since the OH 
emission-line strengths appeared to be much more stable 
during the Mrk 766 observations, leading to a cleaner sky 
subtraction. 

3. STELLAR VELOCITY DISPERSIONS 

To measure stellar velocity dispersions, we used the 
Ca II triplet region (rest frame ~8500A) in the optical 
spectra, and Mg I and the CO line region (rest frame 
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Fig. 1. — Velocity dispersion measurements using the Palomar 
DBSP and the Keck ESI spectra. The region including the main 
stellar features around the Ca II triplet is shown together with the 
best-fit template (thick red line). 

~ 1.5-1. 6 ^tm) in the near-IR spectra. Velocity disper- 
sions were determined by comparing in pixel space the 
observed galaxy spectra with stellar templates convolved 
with a kernel representing Gaussian broadening. Veloc- 
ity template stars of various spectral types were observed 
with the same instrumental setup at each telescope, mini- 
mizing any systematic errors due to the different spectral 
resolution of various observing modes. The minimum \ 2 
fit was performed using the Gauss-Hermite Pixel Fitting 
software (van der Marel 1994). The merit of fitting in 
pixel space as opposed to Fourier space is that typical 
AGN narrow emission lines and any residuals from night- 
sky subtraction can be easily masked out (e.g., Treu et 
al. 2004; Woo et al. 2005; 2006), and the quality of the 
fit can be directly evaluated. We used low-order poly- 
nomials (order 2-4) to model the overall shape of the 
continuum. Extensive and careful comparison was per- 
formed using continua with various polynomials, and the 
measured velocity dispersions based on each continuum 
fit with various order polynomials are consistent within 
the errors (see e.g., Woo et al. 2004; 2005; 2006). Thus, 
continuum fitting with low-order polynomials does not 
significantly affect the velocity dispersion measurements. 

3.1. Optical Data 

For galaxies with optical spectra we mainly used the 
strong Ca II triplet lines (8498, 8542, 8662 A) to mea- 
sure the stellar velocity dispersions. We used all veloc- 
ity template stars with various spectral types (G8-K5 
giants) observed with the same instrumental setups, in 
order to account for possible template mismatches. Fits 
were performed with each template and the mean of these 
measurements was taken as the final velocity dispersion 
measurement. The rms scatter around the mean of in- 
dividual measurements was added in quadrature to the 
mean of the individual measurement errors from the x 2 
fit. Given the high S/N of the observed spectra, the un- 
certainties are typically ~5%. 

We were able to measure the stellar velocity dispersion 
of eight objects: SBS 1116+538A, Arp 151, Mrk 1310, 
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Fig. 2. — Galaxies with contaminated stellar lines. These galax- 
ies clearly show the Ca II triplet lines; however, due to AGN con- 
tamination, the line profiles are asymmetric. Also, the presence of 
various AGN emission features complicates continuum subtraction. 
The spectra for the first three objects (Mrk 871, NGC 4748, Mrk 
766) are from Palomar DBSP data. The data for Mrk 290 data 
were obtained at the Lick telescope, and Mrk 142 data are from 
SDSS. 

Mrk 202, IC 4218, NGC 5548, NGC 6814, and MCG- 
6-30-15. Among these, IC 4218 and MCG-6-30-15 were 
excluded in the A/bh — c* relation analysis since we were 
not able to measure the reverberation time scales (Bcntz 
et al. 2009b). Figure 1 shows examples of the best- fit 
template along with the observed galaxy spectra. Our 
measurements based on high-quality data are generally 
consistent within the errors with previous measurements 
from the literature for these galaxies, although many of 
the previously existing measurements had very large un- 
certainties (see £)3. 31 below). 

In the case of the other five objects (Mrk 871, NGC 
4748, Mrk 766, Mrk 290, and Mrk 142), the Ca II triplet 
lines were contaminated by AGN emission lines as shown 
in Figure 2; hence, we did not obtain reliable stellar ve- 
locity dispersion measurements. Three of these objects, 
that have black hole mass measurements from reverber- 
ation mapping, were observed in the near-IR, and we 
determined stellar velocity dispersions for two of them 
as described in the next section. 

3.2. Infrared Data 

Because of the difficulty of measuring the stellar veloc- 
ity dispersion in the optical, we obtained near-IR spectra 
for three galaxies with reverberation black hole masses 
(NGC 4748, Mrk 766, Mrk 142) using OSIRIS. We also 
observed several velocity template stars with K and M 
spectral types. As noted by other studies (e.g., Wat- 
son et al. 2008), M giant stars give better fits com- 
pared to K giants since M giants are the representative 
stellar population in the near-IR bands. After experi- 
menting with various spectral type velocity templates, 
we decided to use M1III and M2III stars for velocity dis- 
persion measurements and took the mean as the final 
measurement with an uncertainty given by the quadra- 
ture sum of the fitting error and the rms scatter of re- 
sults obtained from different templates. Figure 3 shows 
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Mrk 142 

1.5 1.55 1.6 

Rest wavelength (yitm) 

Fig. 3. — Velocity dispersion measurements using the OSIRIS 
spectra. The region including the main stellar features is shown 
together with the best-fit template (thick red line). Individual stel- 
lar lines are identified in the first panel. The regions around narrow 
AGN emission lines and regions of severe template mismatch, iden- 
tified by shaded bands, are masked out before fitting. Residuals 
from the fit are shown at the bottome of NGC 4748 and Mrk 766. 

the observed galaxy spectra compared with the best-fit 
template. NGC 4748 shows a relatively good fit, while 
Mrk 766 presents stronger residuals due to template mis- 
match. All major if-band stellar lines, such as Mg 1 1.488 
fim, Mg I 1.503 fim, CO(3-0) 1.558 (im, CO(4-l) 1.578 
/xm, Si I 1.589 /im, and CO(5-2) 1.598 ^m, are present 
in both objects. In the case of Mrk 142, which is fainter 
and at higher redshift than the other two objects, stellar 
absorption lines cannot be clearly identified and we could 
not obtain a reliable template fit or velocity dispersion 
measurement. 

3.3. Previous Velocity Dispersion Measurements 

Several objects in the LAMP sample have preexisting 
velocity dispersion measurements in the literature. In 
most cases the published results have fairly large uncer- 
tainties. Here we briefly summarize the previous mea- 
surements and describe the available spectroscopic data 
from the Sloan Digital Sky Survey (SDSS). 

SBS 11 16+ 583 A: IGreene fc Hoi (f2006h measured a = 
50.3 ± 18.0 km s _1 using the SDSS spectrum. However, 
the velocity dispersion was not clearly resolved by the 
SDSS data, which have an instrumental dispersion of ~ 
70 km s" 1 . 

Arp 151: IGreene fe Hoi ((20061 ) measured cr = 124 ± 12 
km s _1 from the SDSS spectrum. 

Mrk 1310: Using the SDSS data. IGreene"fcHol ((20061 ) 
reported a — 50 ± 16 km s _1 , which is below the resolu- 
tion limit of the SDSS spectrum. 

Mrk 202: IGreene fc Hoi (|2006l) measured a = 85.6 ± 14 
km s _1 from the SDSS spectrum. 

Mrk 290: In the SDSS spectrum, the Ca II triplet lines 
are only weakly detected, and they appear to be badly 
contaminated by Ca II emission. The blue region of the 
spectrum is too dominated by the AGN for any stellar 
features to be clearly visible. 



NGC 4253 (Mrk 766): IBotte et all ((20051) reported a 
dispersion of cr = 81 ± 17 from the Ca II triplet lines. 
However, the Ca II triplet region is moderately contami- 
nated by emission lines. Their dispersion measurements 
were carried out by cross-correlating the galaxy spectrum 
with that of a template star. Without a direct fit of the 
broadened template star to the galaxy spectrum, it is not 
clear how badly the measurement may be affected by the 
emission-line contamination. 

An SDSS spectrum of Mrk 766 exists, but it appears 
to have been taken with the SDSS fiber positioned about 
3" away from the galaxy nucleus. Since the SDSS fiber 
missed the galaxy nucleus, the spectrum appears similar 
to that of a Seyfert 2 galaxy, dominated by starlight and 
with very weak broad components to the Balmer lines. 
Our own data from the Lick monitoring campaign, on the 
other hand, show a strongly AGN-dominated spectrum 
with much more prominent broad lines. Thus, while 
the Ca II triplet and other stellar absorption lines are 
clearly detected in the SDSS data, the offset position of 
the SDSS spectrum means that it should not be used to 
infer the bulge velocity dispersion. 

NGC 5548: The most recent measurements of the ve- 
locity di spersion for this gal axy include a = 180 ± 6 
km s" 1 dFerrarese et aill2001[ ) and cr = 201 ± 12 km s" 1 
(|Nelson et al.ll2004D . 

NGC 681 4: INelson fc Whittle (fl995l) listed a velocity 
dispersion of 115 ± 18 km s _1 . 

Mrk 142: The SDSS spectrum of Mrk 142 (SDSS 
J102531. 28+514034. 8) shows that the object is strongly 
AGN dominated, and the stellar Ca II triplet lines are 
severely contaminated by Ca II emission. In addition, the 
spectral region around Mg b and Fe 5270 is completely 
dominated by strong Fe II emission, precluding any pos- 
sibility of measuring the stellar velocity dispersion in the 
optical. 

MCG-6-30-15: IQliva et all (fl999l) measured a = 159 
km s~ x from template fits to the CO and Si features 
in the H band, but they did not list estimates o f the 
measurement uncertainties. Mc Hardv et "all ([2005) used 
a spectrum around the Ca II triplet to measure cr = 
93.5 ± 8.5 km s , using a cross-correlation technique. 

All of the previous measurements with sufficiently 
high instrumental resolution are consistent with our new 
measurements within the errors, except for one object, 
MCG-6-30-15. The previously reported Ca Il-based ve- 
locity dispersion and the -ff-band velocity dispersion of 
MCG-6-30-15 are dramatically different, and our mea- 
sured dispersion based on the Ca II triplet is significantly 
lower than either previous result. Unfortunately, we can- 
not determine the origin of the discrep ancy other than 
noting that the Ca II dispersion from fMcHard v et al.l 
(2005) was determined using a cross-correlation method, 
and therefore it is not easy to directly compare the qual- 
ity of the result with our measurement. As an additional 
check, we carried out an independent measurement of cr» 
from t he Palomar DBSP data using the template-fitting 
code of lBarth et aT1(|2002D . This gave cr* = 75±5 kins' 1 , 
con sistent with the value of 76 ±3 km s _1 measured with 
the Ivan der Marell (fl99l code. Since MCG-6-30-15 is 
not included in the Mbh — relation analysis in the 
next section, this does not affect our main conclusions. 
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TABLE 2 

Black hole masses and stellar, velocity dispersions 



Galaxy 


VP 

K ne iW/G) 
1O 6 M 


VP ref. 


cr * 
km s~ 1 


error 


cr * 


ref. 


log Mbh/M 


LAMP Sample 


Arp 151 


i 99 +U.lt> 
-■-■ zz -0.22 


5 


118 


4 


this 


work 


6.81 ± 0.12 


IC 4218 






93 


4 


this 


work 




MCG-6-30-15 






76 


3 


this 


work 




Mrk 1310 


41 + ' 12 


5 


84 


5 


this 


work 


6.33 ± 0.17 


Mrk 142 


40+ ' 12 


5 










6.32 ± 0.17 


Mrk 202 


26+ - 13 


5 


78 


3 


this 


work 


6.13 ± 0.22 


NGC 4253 (Mrk 766) 


S2+ 21 
u - JZ -0.20 


5 


93 


32 


this 


work 


6.23 ± 0.30 


NGC 4748 


47+°- 16 
u - 4 '-0.21 


5 


105 


13 


this 


work 


6.39 ± 0.20 


NGC 5548 




6 


195 


13 


this 


work 


7.80 ±0.10 


NGC 6814 




5 


95 


3 


this 


work 


7.25 ± 0.12 


SBS 1116+583A 


l.UO_ 2g 


5 


92 


4 


this 


work 


6.74 ±0.16 


Previous Reverberation Sample 


Ark 120 


27-2+35 


1 


221 


17 




1 


8.15 ± 0.11 


3C 120 


io.i±|;I 


1 


162 


20 




2 


7.72 ± 0.23 


3C 390.3 


52+ 11 ' 7 


1 


273 


16 




1 


8.44 ±0.14 


MRK 79 


q co+2.61 
9 - oz -2.61 


1 


130 


12 




1 


7.70 ±0.16 


MRK 110 


4.57+^ 


1 


91 


7 




3 


7.38 ± 0.14 


MRK 279 


6.35+^ 


1 


197 


12 




1 


7.52 ± 0.15 


MRK 590 




1 


189 


6 




1 


7.66 ± 0.12 


MRK 871 


8.98±lt 


1 


120 


15 




1 


7.67 ±0.12 


NGC 3227 


7.6711* 


1 


136 


4 




1 


7.60 ± 0.24 


NGC 3516 


7.76l 2 2 ;i 


1 


181 


5 




1 


7.61 ± 0.18 


NGC 3783 


r 42 +0.99 
°-* z -0.99 


1 


95 


10 




4 


7.45 ± 0.13 


NGC 4051 


9O7+0.09 
•^°'-0.12 


2 


89 


3 




1 


6.18 ± 0.19 


NGC 4151 


o O1+1-04 

°- di -0.85 


3 


97 


3 




1 


7.64 ±0.11 


NGC 4593 


i 7O+0.38 


4 


135 


6 




1 


6.97 ±0.14 


NGC 7469 


9 91 +0.25 
z - zl -0.25 


1 


131 


5 




1 


7.06 ± 0.11 


PC 1426+215 


2361™ 


1 


217 


15 




5 


9.09 ± 0.16 



Note. — Col. (1) name. Col. (2) virial product (M Bi i=fx VP) based on the line dispersion 
("line) from reverberation mapping. Col. (3) reference for virial product. 1. Peterson et al. 2004; 
2. Denney ct al. 2009; 3. Bcntz ct al. 2006; 4. Denney et al. 2006; 5. Bentz et al. 2009b; 6. 
weighted mean of Bcntz et al. 2007 and Bcntz ct al. 2009b. Col. (4) stellar velocity dispersion. 
Col. (5) error of stellar velocity dispersion. Col. (6) reference for stellar velocity dispersion. 
1. Nelson et al. 2004; 2. Nelson & Whittle 1995; 3. Ferrarese et al. 2001; 4. Onken et al. 
2004; 5. Watson et al. 2008. Col. (7) black hole mass calculated assuming the virial coefficient, 
log / = 0.72 ±0.10. The uncertainty in the black hole mass is calculated by adding in quadrature 
the measurement uncertainty of the virial product in logarithmic space and the uncertainty in the 
virial coefficient (0.1 dex) as e log Mbh = \/ eVP 2 /(VP lnlO) 2 + 0.1 2 . The average of the positive 
and negative errors on the virial product is taken as the symmetric error (eVP). 



4. THE M BH - cr. RELATION OF ACTIVE GALAXIES 
The Mbh — cr* relation is expressed as 

log(M BH /M Q ) = a + /31og(<7*/200 km s" 1 ). (2) 

For a sample of quiescent galaxies for which the sphere 
of influence of the black hole is spatially resolved in stel- 
lar or gas kinematics studies, Ferrarese & Ford (2005) 
measured a = 8.22 ± 0.06 and f3 = 4.86 ± 0.43 with no 
evidence for intrinsic scatter. Using a larger sample, in- 
cluding late-type galaxies, Giiltekin et al. (2009) recently 
reported a = 8.12 ± 0.08 and (3 = 4.24 ± 0.41 with an 
intrinsic scatter <Ti nt = 0.44 ± 0.06. 

In this section, we investigate the Mbh — cr* relation 
of active galaxies, by combining our 8 new velocity mea- 
surements of the LAMP sample with 16 existing mea- 
surements of reverberation black hole mass and velocity 
dispersion taken from the literature. The relevant prop- 



erties of the resulting sample of 24 AGNs are listed in 
Table 2, along with the original references for the mea- 
surements. Note that the selection bias investigated by 
Lauer et al. (2007) is not relevant to these reverberation 
objects since the bias mainly affects black holes in high- 
mass galaxies where the host galaxy luminosity function 
is steeply falling. 

As a first step we determine the slope of the Mbh — cr* 
relation (§4.1) assuming that the virial coefficient / is 
unknown but does not depend on black hole mass. We 
also determine the intrinsic scatter in the relation. Then 
in § 4.2 we determine the mean value and intrinsic scatter 
of the virial coefficient / by assuming that the slope and 
zero point of the Mbh — cr* relation are the same for 
active and quiescent galaxies. Note that in the following 
analysis the virial coefficient / and the zero point a are 
degenerate; therefore, additional information is needed 
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to determine what fraction of the intrinsic scatter is due 
to / and how much is due to a, as we discuss below. 

4.1. The Slope 

Active galaxies appear to obey an A/bh — o~* relation 
with a slope consistent with that of quiescent galaxies 
(e.g., Gebhardt et al. 2000; Ferrarese et al. 2001). How- 
ever, the slope has not been well determined so far, 
because of limitations in sample size and in the dy- 
namic range of -Mbh for active galaxies with reverber- 
ation and stellar velocity dispersion measurements. The 
local AGN samples with single-epoch black hole masses 
selected from the SDSS are much larger than the rever- 
beration sample (e.g., Greene & Ho 2006; Woo et al. 
2008), but these mass estimates have larger uncertain- 
ties (~ 0.5 dex) and may suffer from unknown systematic 
errors (Collin et al. 2006). 

By adding our LAMP galaxies to the previous rever- 
beration sample, we significantly increase the sample 
size, particularly at low mass scales (M BH < 10 7 M Q ), 
and extend the dynamic range to almost 3 decades in 
Mbh- This enables substantial progress in determining 
the slope of the Mbh — o~* relation, independently for 
active galaxies. One important caveat is that the virial 
coefficient / may vary across the sample. Unfortunately, 
with existing data we cannot constrain the virial coeffi- 
cient of each object. Therefore, in the analysis presented 
here we will assume that / is independent of black hole 
mass, although we allow for intrinsic scatter at fixed 
black hole mass. In the next section we will take the 
opposite viewpoint and assume that the slope is known 
from quiescent samples when we investigate the intrin- 
sic scatter in / and its dependence on Mbh and AGN 
properties. 

In practice, our observable is the so-called virial prod- 
uct VP= of inc -RBLR/G, which is assumed to be related 
to the black hole mass by a constant virial coefficient 
A^bh^/xVP. Substituting into Equation^ one obtains 

logVP^ + /3 log(^- T )-lo g/ . (3) 

In order to determine the slope /3 of the relation, we 
compute the posterior distribution function P. Assum- 
ing uniform priors on the parameters a' — a — log /, /3, 
and intrinsic scatter eri n t, the posterior is equal to the 
likelihood given by 

- 1 b fs fNvp.-' | -WM' +lt ,^r 

(4) 

where 

e tot,i = e logVP; + e \oga,i + °fnt7 (5) 

and ei g vp is the error on the logarithm of the virial 
product (= eyp/(VP lnlO)) using the average of the pos- 
itive and negative errors on the virial product as eyp, 
^iog (7* is the error on the logarithm of the stellar veloc- 
ity dispersion (= e rT ,/(o'* lnlO)), and a- ul t is the intrinsic 
scatter, which is treated as a free parameter. Note that 
our adopted likelihood is equivalent to that adopted as 
default by Giiltekin et al. (2009) for their maximum like- 



lihood analysis (see their Appendix A. 2), with the in- 
clusion of an additional Gaussian term to describe the 
uncertainties in velocity dispersion. After marginalizing 
over the other variables, the one-dimensional posterior 
distribution functions yield the following estimates of the 
free parameters: a' = 7.28 ± 0.14, /3 = 3.55 ± 0.60, and 
(Tint = 0.43 ± 0.08. 

To facilitate comparison with earlier studies (e.g., 
Tremaine et al. 2002), we also determine the slope with 
a "normalized" x 2 estimator defined as 

i—l £ tot,z 

where the intrinsic scatter <7j n t = 0.43 is set to yield % 2 
per degree of freedom equal to unity. The best estimate 
of the slope based on this normalized \ 2 estimator is 
f3 = 3.72lg3g. If we assume no intrinsic scatter, then the 

best-fit slope becomes slightly steeper as f3 — 3.97lg 
However, confirming our previous result based on a full 
posterior analysis, the intrinsic scatter cannot be negli- 
gible since the \ 2 P er degree of freedom would be much 
larger than unity (~14) if the intrinsic scatter is assumed 
to be zero. 

These results indicate that the slopes of the Mbh — <?* 
relation for active and quiescent galaxies are consistent 
within the uncertainties. The total intrinsic scatter is 
0.43 ± 0.08 dex, which is a combination of the scatter of 
the AIbh — °~* relation and that of the virial coefficient. 
Without additional information, we cannot separate the 
contribution of the intrinsic scatter on the virial coeffi- 
cient (i.e. the diversity of broad line region kinematic 
structure amongst AGN) from the intrinsic scatter in 
black hole mass at fixed velocity dispersion (i.e., the di- 
versity of black hole masses amongst galaxies). We will 
address this issue in §4.2. 

4.2. The Virial Coefficient 

In order to transform the virial product into an actual 
black hole mass, we need to know the value of the virial 
coefficient /. Our current understanding of the geome- 
try and kinematics of the BLR is insufficient to compute 
/ from first principles. This creates a systematic un- 
certainty in the determination of black hole mass from 
reverberation mapping. As an example, if a spherical 
isotropic velocity distribution is assumed, then the veloc- 
ity dispersion of the gas is a factor of larger than the 
measured line-of-sight velocity dispersion (cri ine ) , hence a 
value of / = 3 is often assumed for a spherically sym met- 
ric BLR (e.g.. IWandel et al.|[l999t IKaspi et alboOOfl . If, 
instead, the BLR is described by a circular rotating disk, 
then the value of / can be a few times larger than the 
isotropic case, depending on the inclination angle and the 
scale height of the disk (see Collin et al. 2006 for more 
details). 

As mentioned above, determining the virial coefficient 
for individual AGNs is not feasible due to the limited 
spatial resolution of the current and foreseeable future 
instruments since the angular size of the BLR is mi- 
croarcseconds. Therefore, in practice an average virial 
coefficient is usually determined by assuming that local 
Seyfert galaxies and quasar host galaxies follow the same 
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Mbh — c* relation (Onken et al. 2004) and deriving the 
appropriate value for /. Using this method, Onken et 
al. (2004) measured (/) = 5.5 ± 1.8 based on a sample of 
14 reverberation-mapped Seyfert galaxies with measured 
stellar velocity dispersions. 

In this section we apply the same methodology to our 
enlarged sample to determine the average value of / as 
well as the intrinsic scatter, assuming the intercept (a) 
and slope (/3) of the Mbh — cr* relation for quiescent 
galaxies. In practice, we modify the likelihood given in 
Equation |4] by fixing the slope and intercept to the val- 
ues obtained by two independent groups, Ferrarese & 
Ford (2005) and Giiltekin et al. (2009). Figures 4 and 
5 show the best-fit Mbh — f* relations for the enlarged 
reverberation sample, respectively assuming the quies- 
cent Mbh — v* relation from Ferrarese & Ford (2005) 
and that from Giiltekin et al. (2009). Adopting the 
M B h - a* relation of Ferrarese & Ford (2005; a = 8.22 
and /3 = 4.86), we determine log/ = 0.71 ± 0.10 and 
Omt = 0.46io'o9 dex. Adopting the A/bh — c* relation 
from Giiltekin et al (2009; a = 8.12 and /3 = 4.24), we 
obtain log/ = 0.72±g;5o and cr int = 0.44 ± 0.07. The 
average values are in good agreement with / = 5.5 ± 1.8 
(i.e., log / = 0.74l°;i°) as found by Onken et al. (2004). 
The average value is inconsistent with that expected for 
a spherical BLR (/ = 3), and closer to that expected for 
a disk-like geometry, although there may be a diversity 
of geometries and large-scale kinematics (e.g., Bentz et 
al. 2009b; Denney et al. 2009). 

As discussed at the end of § 4.1, the intrinsic scatter de- 
termined in this section is a combination of intrinsic scat- 
ter in / and in zero point of the Mbh — &* relation (which 
are degenerate in Equation 4). We need additional infor- 
mation to break this degeneracy. If the intrinsic scatter 
of the active galaxy Mbh — &* relation were close to 0.31 
dex (as determined for elliptical galaxies by Giiltekin et 
al. 2009), then the residual intrinsic scatter of the virial 
coefficient would be 0.31 dex. In contrast, if the scatter 
of the Mbh ^ relation of active galaxies were closer 
to 0.44 dex (as found for all galaxies by Giiltekin et al. 
2009), then the intrinsic scatter in / would be close to 
zero. Although we cannot break this degeneracy with 
current data, the bottom line is the same as far as de- 
termining black hole masses from reverberation mapping 
data is concerned (and from single-epoch data, as a con- 
sequence). Since for the purpose of the calibration of / 
the two sources of scatter are indistinguishable, the un- 
certainty in our calibration of / (whether it is due to 
diversity in BLR physics or diversity in the galaxy- AGN 
connection) is the largest remaining uncertainty in black 
hole mass determinations that rely on this technique. 

To investigate any dependence of the virial coefficient 
on AGN properties, we collected from the literature the 
H/3 line widths (Vfwhm and an nc ) and the optical AGN 
luminosities for reverberation-mapped AGNs, corrected 
for the host-galaxy contribution (Peterson et al. 2004; 
Bentz et al. 2009a). When multiple luminosity mea- 
surements are available for given objects, we used the 
weighted mean of the luminosity (Table 9 of Bentz et al. 
2009a). For the seven LAMP objects, with Mbh and ve- 
locity dispersion measurements, Hubble Space Telescope 
imaging is not yet available for measuring the starlight 
correction to the spectroscopic luminosity. Thus we fil- 
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Fig. 4. — Top: Mbh — °"* relation of the combined reverbera- 
tion sample, adopting the slope from Ferrarese & Ford (2005) and 
the best-fit virial coefficient from this work. Bottom: Posterior 
probability contour levels (68% and 95%) for the virial coefficient 
/ and its intrinsic scatter, as determined assuming the slope from 
Ferrarese & Ford (2005). The horizontal lines represent the previ- 
ous estimate from Onken et al. (2004; dashed line) with its 68% 
uncertainties (dotted lines). 

ferred the AGN luminosity from the measured time lag 
using the size-luminosity relation (Bentz et al. 2009a). 
The bolometric luminosity was calculated by multiplying 
the optical AGN luminosity by a factor of ten (Woo & 
Urry 2002), and we assumed 0.3 dex error in the Edding- 
ton ratio. 

In Figure 6, we plot the residuals from the Mbh — o* 
relation with a = 8.12, /3 = 4.24 taken from Giiltekin 
et al. (2009), and log/ = 0.72 as determined in §6.2. 
The residuals were determined as A log Mbh = log VP + 
log /— (a+(3 log <7* ) . To test for potential correlations be- 
tween residuals and line properties or Eddington ratios, 
we fit each data set with a straight line by minimizing 
the normalized \ 2 estimator. We find no dependence of 
the residuals on the line width (either Vfwhm or crime). 
In the case of the Eddington ratio, we find a weak de- 
pendence with a non-zero slope ( — 2.3^5 l)' significant at 
the 95% level; however this weak trend is mainly due to 
one object with the lowest Eddington ratio in the sample. 
When we remove that object, the best fit slope becomes 
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Fig. 5. — As in figure [4] but adopting the slope from Giiltekin et 
al. (2009). 

consistent with no correlation. Further investigation us- 
ing a larger sample, evenly distributed over the range of 
the Eddington ratio, is required to better constrain the 
residual dependence on the Eddington ratio. In the case 
of the line width, it is not straightforward to test the 
dependence since the line width of most of our objects 
is relatively small, crime < 2000 km s _1 . By dividing our 
sample at erii nc = 1500 km s" 1 into two groups of similar 
sample size, we separately measured the virial coefficient 
for narrower-line and broader-line AGNs. The difference 
in the virial coefficient is A log/ = 0.1 — 0.2, which is 
not significant given the uncertainty of 0.15 dex on the 
virial coefficient. 

5. DISCUSSION 

We present the Mbh — cr* relation of the reverberation 
sample in Figure 7, using the slope, /3 = 3.55 ± 0.60, 
determined in §4.1 and the average virial coefficient, 
(log/) = 0.72 ±0.10, determined in §4.2. Compared 
to the local quiescent galaxies, active galaxies follow a 
consistent Mbh — cr* relation with a similar slope and 
scatter. Note that the mean black hole mass of the re- 
verberation sample ((logMBH/M Q ) = 7.3 ± 0.74) is an 
order of magnitude smaller than that of quiescent galax- 
ies ((log Mbh/Mq) = 8.2 ±0.79). The slightly shallower 
slope of the reverberation sample is consistent with the 
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Fig. 6. — Dependence of residuals from the Mbh — cr* relation 
(A log M BH = (log VP + log/) - (a + /31ogcr/(200 km s" 1 )) on 
parameters related to the accretion state: VpwHM/cime (top left); 
Eddington ratio (top right); Vfwhm (bottom left); line dispersion 
c line (bottom right) of the H/3 line. In this plot we adopt the local 
relation with a = 8.12, = 4.24 taken from Giiltekin et al. (2009), 
and log / = 0.72 as determined in §6.2. 

trend in the quiescent galaxies that the slope is shallower 
for galaxies with lower velocity dispersion (cr* < 200 
km s^ 1 ); see Giiltekin et al. (2009). In contrast, the slope 
of late-type quiescent galaxies (4.58 ±1.58) seems higher 
than the slope of our active galaxies, which are mainly 
late-type galaxies. However, given the uncertainty in the 
slope (P = 3.55±0.60) of the active galaxy Mbh — cr* rela- 
tion, the difference between quiescent and active galaxies 
is only marginal. We did not attempt to divide our sam- 
ple into various morphology groups or a few mass bins 
to test the dependence of the slope, since the sample size 
is still small and biased toward lower mass objects. 

The smaller mean Mbh of the reverberation sample is 
due to the relative difficulty of measuring the stellar ve- 
locity dispersions for the more luminous and higher red- 
shift quasars that have the largest black hole masses. Al- 
though there are 17 quasars with measured reverberation 
black hole masses, the stellar velocity dispersions have 
not been measured for most of those high-luminosity ob- 
jects (see Watson et al. 2008; Dasyra et al. 2007). In our 
current sample, there are only three objects above 3 x 10 s 
Mq . In order to better constrain the Mbh — cr* relation 
of active galaxies, it will be necessary to pursue further 
measurements of velocity dispersions for reverberation- 
mapped AGNs at high masses. Additional stellar veloc- 
ity dispersion measurements will be available based on 
Keck OSIRIS spectra in the near future (Woo & Malkan, 
in preparation). 

The key assumption in determining the slope of the 
Mbh — a* relation is that the virial coefficient does not 
systematically vary as a function of black hole mass. 
Currently, we do not have sufficient information to test 
whether this assumption is valid. In an empirical study 
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Fig. 7. — The Mbh — °* relation of active galaxies with reverberation black hole masses (blue), compared with non-active galaxies with 
dynamical black hole masses. The reverberation masses were determined assuming the virial coefficient, log / = 0.72 ± 0.10. The solid line 
is the best-fit slope of the active galaxies while the dashed line is the best fit to the inactive galaxy samples from Giiltekin et al. (2009). 



of the systematic effects of emission-line profiles, Collin 
et al. (2006) showed that it was necessary to use a 
smaller value of the virial coefficient for AGNs with larger 
Vfwhm- In contrast, the change of the virial coefficient 
was not required if the line dispersion (ay lae ) was used 
instead. These results reflect the fact that there is a 
large range of Vfwhm /ciine ratio, and that systematic 
uncertainties are higher if Vfwhm is used. 

With the enlarged sample of 24 AGNs, we investigated 
the dependence of the virial coefficient on the line proper- 
ties and the Eddington ratio. However, we did not clearly 
detect any increasing trend of the intrinsic scatter with 
the line profile, line width, or Eddington ratio, implying 
that the Mbh based on the mean virial coefficient is not 
systematically over- or under-estimated as a function of 
line width or the Eddington ratio. We note that since 
our sample is still small and biased to narrower-line and 
lower Eddington ratio objects, further investigation with 
a larger sample, particularly at higher mass scales, is re- 



quired to better constrain the systematic errors. 

For single-epoch black hole mass determination, we 
suggest using log/ = 0.72 ± 0.10 when the gas veloc- 
ity is determined from the line dispersion (<7ii nc ). If the 
FWHM of an emission line (Vfwhm) is measured instead 
of the line dispersion, then additional information on the 
Vfwhm /fiine ratio is needed in the black hole mass cal- 
culation. The value of VFWHM/cii no is typically assumed 
to be 2 (Netzer 1990), although there is a range of values. 
For example, the mean of VpwHM/^line of our reverber- 
ation sample is 2.09 with a standard deviation of 0.45 
(see also Peterson et al. 2004 and McGill et al. 2008). 
We did not attempt to determine an average virial co- 
efficient using the virial product based on Vfwhm- A 
correction factor depending on VFWHM/Viine should be 
applied in addition to the virial coefficient, log / = 0.72, 
when gas velocity is determined from Vfwhm- Our em- 
pirical measurement of the mean virial coefficient is not 
very different from that of Onken et al. (2004) although 
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the sample size and the dynamical range of the reverber- 
ation sample is significantly improved. 

We note that our results on the mean virial coefficient 
are based on broad-line dispersions measured from the 
rms spectra, which are not generally available for single- 
epoch black hole mass determinations. Detailed com- 
parison of broad-line profiles and widths between rms 
and single-epoch spectra is necessary to understand the 
additional uncertainty of single-epoch black hole mass 
estimates (see Denney et al. 2009), and we will pursue 
further investigations of this issue using LAMP data in 
the future. 

6. SUMMARY 
We summarize the main results as follows. 

(1) Using the high S/N optical and near-IR spectra 
obtained at the Keck, Palomar, and Lick Observatories, 
we measured the stellar velocity dispersion of 10 local 
Seyfert galaxies, including 7 objects (SBS 1116+ 583A, 
Arp 151, Mrk 1310, Mrk 202, NGC 4253, NGC 4748, 
NGC 6814) with newly determined reverberation black 
hole masses from the LAMP project (and NGC 5548, IC 
4218 and MCG-6-30-15). 

(2) For a total sample of 24 local AGNs, combin- 
ing our new stellar velocity dispersion measurements of 
the LAMP sample and the previous reverberation sam- 
ple with the stellar velocity dispersion from the liter- 
ature, we determined the slope and the intrinsic scat- 
ter of the A/bh — cr* relation in the range of black hole 
mass 10 6 < Mbh/M q < 10 9 . The best- fit slope is 
P = 3.55 ± 0.60, consistent within the uncertainty with 
the slope of quiescent galaxies. The intrinsic scatter, 
0.43 ± 0.08 dex, of active galaxies is also similar to that 
of quiescent galaxies. Thus, we find no evidence for de- 
pendence of the present-day Mbh ~ c* relation slope on 
the level of activity of the central black hole. 

(3) We determined the virial coefficient using the slope 
and the intercept of the Mbh — &* relation of quies- 
cent galaxies, taken from two groups (Ferrarese & Ford 
2005; Giiltekin et al. 2009). The best-fit virial coefficient 
(log/) is 0.71 ± 0.10 (0.72l°;° 9 ) with an intrinsic scat- 
ter, 0.46^09 (0-44 ± 0.07) with the slope, 4.86 (4.24) 
taken from Ferrarese & Ford (Giiltekin et al. 2009). We 



take / = 5.2 (i.e., log / = 0.72) as the best value of the 
mean virial coefficient, which is a factor of ~1.7 larger 
than the standard factor obtained for isotropic spatial 
and velocity distribution. The substantial intrinsic scat- 
ter indicates that the virial coefficient is the main source 
of uncertainties in determining black hole masses, either 
using reverberation mapping or single-epoch spectra. 
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